Husimi functions at dielectric interfaces: Inside-outside duality for optical systems and beyond 
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We introduce generalized Husimi functions at the interfaces of dielectric systems. Four different functions 
can be defined, corresponding to the incident and departing wave on both sides of the interface. These functions 
allow to identify mechanisms of wave confinement and escape directions in optical microresonators, and give 
insight into the structure of resonance wave functions. Off resonance, where systematic interference can be 
neglected, the Husimi functions are related by Snell's law and Fresnel's coefficients. 

PACS numbers: 03.65.Sq, 03.65.Nk, 42.25.-p, 42.60.Da 
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Optical microresonators receive growing interest over the 
last years, because of the intricate interplay of shape (leading 
to irregular classical ray dynamics), openness of the system 
(offering means of excitation and escape), and the wave na- 
ture of the field. This interplay, together with the promising 
prospect of applications in future communication devices, has 
stimulated experiments [|lj [| ||] as well as theoretical investi- 
gations [Q, ||] that were based on concepts well-known from 
scattering theory, classical ray dynamics, semiclassics, and 
quantum chaos [||] . A particularly useful tool to study waves 
in dynamical systems is the Husimi representation of the wave 
function in classical phase space. So far, the Husimi represen- 
tation was mostly used to study the closed analogues of opti- 
cal microsystems, with the dielectric interface being replaced 
by hard walls, and the principal confinement and radiation di- 
rections have been inferred by adding the laws of reflection 
and refraction by hand. Most notably, efforts in this direction 
suffer from the fact that the incident and emerging wave com- 
ponents cannot be discriminated by the conventional Husimi 
representation. The reasons for these fundamental shortcom- 
ings arise from the facts that the wave function of the dielec- 
tric system is only partially confined by the internal reflection 
at the refractive index boundary, and that it is affected by the 
different nature of the boundary conditions, which are neither 
of Dirichlet nor of von-Neumann type but of a mixed type 
that follows from Maxwell's equations (both the wave func- 
tion and its derivative are non- vanishing at the interface). 

In this paper we introduce four Husimi functions appropri- 
ate for dielectric interfaces, corresponding to the intensity of 
incident and emerging waves at both sides of the interface. 
In the regime of ray optics it will be demonstrated that these 
Husimi functions are related across the interface via Fres- 
nel's formulas, with phase space being deformed according 
to Snell's law. This connection can be seen as a new vari- 
ant of the inside-outside duality . However, ray optics only 
applies when the wave length is short and when systematic in- 
terference effects can be neglected. The Husimi functions do 
not require these limitations and develop their full predictive 
power especially when systematic interference effects lead to 
strong deviations from Fresnel's and Snell's laws. In partic- 




FIG. 1 : Refractive-index boundary of a dielectric system. 



ular, ray optics breaks down close to resonances, where the 
internal part of the scattering wave function is known to be 
almost independent of the incoming wave that excites the sys- 
tem. The Husimi functions still provide an accurate represen- 
tation of the wave function (in particular, they nicely display 
the radiation directions of the field). We illustrate these fea- 
tures using the dielectric circular disk and an annular system 
as examples H . 

Microresonators and ray optics. Consider the light that 
illuminates and permeates a dielectric system confined by a 
refractive-index boundary, as shown in Fig. [j] The disk (and 
also the annulus) are confined by a circular dielectric inter- 
face of radius r c = 1. Polar coordinates r, will be used to 
parameterize position space. We distinguish four wave com- 
ponents: The incident (inc) wave and the emerging (em) wave 
on both sides of the interface (region with refractive index 
no = 1 outside the microresonator, region 1 with n\ = n 
inside the microresonator). In the regime of ray optics, the 
wave is represented by rays, where the angles of incident and 
emerging rays are related by the laws of reflection and Snell's 
law, n sin xi = sin \o- For our circular interfaces Snell's law 
is equivalent to conservation of the angular-momentum vari- 
able m = kj sin Xj , where ko and k\ = nko are the wave 
number in each region. The ray intensities on either side fol- 
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The reflection and transmission coefficients Ri and Ti are re- 
lated by 

Ro(xo) = Ri(xi) = R , Mxo) = Ti(xi) = 1-R, 

(3) 

with [§ ii = sin 2 (xi - Xo)/sin 2 (xi + Xo). 

Husimi functions at a dielectric interface. The Husimi 
function at the system boundary of closed systems was in- 
troduced in Refs. [^[ [l(J by projection of the conventional 
Husimi function from full phase space [coordinates (r, <£), 
momentum (kj sinxj, kj cosxj)] onto the reduced phase 
space at the boundary r = r c with coordinates (j) = $ and 
sin Xj JTl| ] . The four different Husimi functions (correspond- 
ing to the incident and emerging wave at both sides of a di- 
electric interface) can be constructed by the same procedure 
when the appropriate boundary conditions are employed. The 
intensities 7(0, sin x) = H (0, sin x) d<p d sin x will turn out 
to be related by the laws (Jfy when ray optics applies, but 
also accurately describe the wave function when ray-optical 
relations across the interface break down, as expected, e.g., 
for resonances. 

The conventional Husimi function for a given wave func- 
tion ^(r, $) of the dielectric system is obtained as the overlap 
with a wave packet with minimal uncertainty in the variables 
(r, <E>) for real space and (kj sin Xj , kj cos Xj) for momentum 
space. The projection onto the boundary can be formulated 
rigorously fllCQ: The wave function is expressed by means 
of advanced and retarded Green's functions, which in turn 
allow to distinguish between incident and emerging waves. 
Green's formula is then used to express the solution tyj of the 
Helmholtz equation in region j = 0, 1 as an integral over the 
boundary, involving both \£ j and its normal (radial) derivative 
fyj . A semiclassical (saddle-point) approximation then allows 
to identify in these expressions the following four different 
Husimi functions on the interface, 
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with the angular-momentum dependent weighting factor 
Tj — ^Jrij cos Xj Jl2| ] • Here the functions 
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are overlaps with the minimal-uncertainty wave packet 

£(r,$;0,sinx) = ^ e -^(Mi-^^ksin X ( m i) 



x(air) 4 S(r — r c ) 



(a periodic function in <1>), which is restricted to the interface 
and centered around ((/>, sinx). The parameter a controls its 
extension in ^-direction, thereby also fixing the uncertainty in 
sinx- We set a = \f2jk\. The scaling with k\ results in the 
same resolution in (j) for all four Husimi functions. 

Inside -outside duality. As a consequence of the bound- 
ary conditions derived from Maxwell's equations we find the 
identities fto(0,sinxo) = M^ sm Xi)> K^^^Xo) = 
h[((j), sinxi), where the angles Xi are related by Snell's law. 
From these relations it follows that our Husimi functions 
strictly fulfill the condition of intensity conservation, 



ntf em ((/>,sinxo) +ir 1 em (0,sin X i) 
n^ c ((/>,sinxo)+^r(^smxi) 



(8) 



where the factor n = d sin Xo/d sin xi accounts for the phase- 
space deformation by Snell's law. Additional relations be- 
tween the Husimi functions can be anticipated in the regime of 
ray optics: The intensities on one side of the interface should 
be related to the intensities on the other side by Eqs. (jl|), (||). 
We then expect validity of the resulting inside-outside duality 
relations 
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which express the Husimi functions in region by the Husimi 
functions in region 1. The notation S(-) indicates that the 
approximation is of semiclassical (short wave length) nature; 
most noticeably, the Husimi functions (intensities) are added 
incoherently. The duality relations can also be inverted, 
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However, the Husimi functions in region 1 can only be recon- 
structed in the strip | sin xi | < 1/n, because the rest of phase 
space is isolated from region by total internal reflection. 

The duality relations are exactly fulfilled in two simple 
cases, namely, if one incident or emerging wave vanishes or 
when the two incident waves have the same intensity (the two 
emerging waves then have the same intensity, as well). We 
now test the duality relations in more general situations. 

The panels on the left in Fig. || show the exact Husimi func- 
tions from Eq. (^) for the case of the circular disk which is ex- 
cited by a plane wave at an off-resonant excitation frequency. 
The illuminating plane wave is clearly visible in the Husimi 
function H™, around the polar angle (j) = fan = tt, while the 
focal point of the dielectric disk results in a bright spot in Hq™ 
that is located around <p = 0. There is a close correspondence 
between the Husimi functions of the incident and emerging 
waves, and the deformation of phase space by the stretching 
factor n of Snell's law is clearly visible. 
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FIG. 2: Exact (left) and reconstructed (right) Husimi functions for 
the circular dielectric disk (n = 2) illuminated by a plane wave at 
an off-resonant frequency. The top panel shows the scattering wave 
function in real space. The exact Husimi functions are obtained 
from Eq. (W). The reconstructed Husimi functions are obtained by 
Eqs. (P)-(|l^). Negative Husimi densities are shown in black. The 
dashed lines in the panels for region 1 mark the critical angle of in- 
cidence for total internal reflection. 



The right panels of Fig. || show for comparison the pre- 
dictions of Eqs. (ph-(12). In the reconstruction we used the 
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Same as Fig. H, but for illumination at a resonance fre- 



slightly modified semiclassical versions of Fresnel's coeffi- 
cients and Snell's laws given in Ref. [Jl3|], which are appro- 
priate for the present case of a curved interface (this results in 
a slight, but still noticeable quantitative improvement of the 
reconstruction). We observe a good qualitative and quantita- 
tive agreement with the exact Husimi functions. Regions with 
unphysical negative intensities are small. The most interest- 
ing deviations between the exact and the reconstructed Husimi 
functions occur around the central spot at <p = 7r, sin Xj < = 0, 
where the incoherent predictions of Eqs. (p|)-(|l2|) underesti- 
mate the exact Husimi densities Hq uc , Hq™, while they over- 
estimate the intensities H™ , Hf m in the same area of phase 
space. These deviations arise from a Fabry-Perot like system- 
atic interference which is constructive in backward direction 
at the presently chosen frequency. At other frequencies the in- 
terference is destructive, and both cases alternate periodically. 

Resonances. Figure [^displays the situation for excitation at 
a frequency which is close to a narrow resonance, a whisper- 
ing gallery (WG) mode located around sinxi = 0.667. The 
top panel shows that the wave function is now well confined 
inside the disk (region 1). Correspondingly, the Husimi func- 
tions H™ c and Hf m noticeably exceed the Husimi functions 
Hq 10 and i^o m - Moreover, the Husimi functions H™ c and 



jjem are dominated by me characteristics of the resonance 
wave function and consequentially are almost independent of 
the choice of the exciting wave. (The remnants of the exciting 
plane can be identified when comparing Fig. [| with Fig. ||.) 
Hence the reconstructed Husimi functions deviate noticeably 
from the exact Husimi functions around | sin xi\ ~ 1/n. This 
is no surprise since resonances are formed by systematic con- 
structive interference, and incoherent ray optics cannot be 
expected to apply under these circumstances. Most impor- 
tantly, by principle, the confined wave intensity in the region 
| sin xi | > V n cannot be reconstructed because classically no 
refracted rays ever reach this region (which is dark off reso- 
nance). On the other hand, the exact Husimi functions display 
nicely all the features of the resonance wave function in phase 
space. 

Finally, let us illustrate the usefulness of the Husimi func- 
tions (|[) also for a more complex system than the circular disk, 
the annular system formed by regions of different refractive 
indices that are confined by two eccentric disks. The ray op- 
tics in this system corresponds to nonintegrable dynamics in 
phase space, which allows for a multifaceted set of resonance 
wave functions [||. Off resonance (Fig. the scattering wave 
function enters the dielectric system only barely, and the sit- 
uation is similar to the circular disk because the internal disk 
is not explored extensively. At resonance the situation is very 
different. Figure || shows a typical resonance wave function in 
real space and its Husimi representation in phase space. The 
intensity of the resonance wave function is concentrated on 
straight segments which can be identified as a short stable pe- 
riodic trajectory of the corresponding classical ray dynamics. 
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FIG. 4: Same as Fig. g, but for a dielectric annulus (refractive 
indices no = 1 outside, m = 3 in the annulus, n<i — 6 in the 
inner disk, radii r c = l,r' c = 0.6, displaced by S = 0.22). The 
system is illuminated by a plane wave at an off-resonant frequency 
(illumination direction 0m = ty). 



The Husimi functions display a strong intensity exactly in the 
vicinity of this trajectory in classical phase space. 

In conclusion, we introduced four Husimi representations 
of the scattering wave function at the interfaces of dielectric 
microresonators, corresponding to the incident and emerging 
waves at both sides of the interface. These Husimi func- 
tions are easily computed from the wave function and have 
many desirable properties: They are related by the laws of 
Fresnel and Snell in the ray-optics regime (i.e., short and off- 
resonant wavelength) and generally provide valuable detailed 
insight into the wave dynamics in complex dielectric systems, 
most notably even close to resonances where ray optics breaks 
down. 

We gratefully acknowledge helpful discussions with Jan 
Wiersig, Arnd Backer, and Christian Miniatura. 
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Same as Fig. W, but for illumination at a resonance fre- 
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